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ABSTRACT 

Spectroscopic long-slit observations of the dwarf Irr galaxy IC 10 were conducted at the 6-m 
Special Astrophysical Observatory telescope with the SCORPIO focal reducer. The ionized-gas 
emission spectra in the regions of intense current star formation were obtained for a large number 
of regions in IC 10. The relative abundances of oxygen, N + , and S + in about twenty HII regions 
and in the synchrotron superbubble were estimated. We found that the galaxy-averaged oxygen 
abundance is 12 + log(0/H) = 8.17 ± 0.35 and the metallicity is Z = 0.18 ± 0.14Z Q . Our 
abundances estimated from the strong emission lines are found to be more reliable than those 
obtained by comparing diagnostic diagrams with photoionization models. 

Subject headings: Irr galaxies, IC 10, interstellar medium, chemical composition. 



INTRODUCTION 

The dwarf Irr galaxy IC 10 is a unique testing 
ground for investigating the interstellar medium 
in regions of violent star formation. The galaxy's 
images in the Ha and [SII] lines appear as a single 
giant complex of multiple shells and supershells 
with sizes from 50 to 800-1000 pc (Zucker 2000; 
Wilcots and Miller 1998; Gil de Paz et al. 2003; 
Leroy et al. 2006; Chyzy et al. 2003; Lozinskaya 
et al. 2008). Its stellar population suggests a re- 
cent starburst (t=4 - 10 Myr) and an older star- 
burst (t > 350 Myr) (Hunter 2001; Zucker 2002; 
Massey et al. 2007; Vacca et al. 2007; and ref- 
erences therein). The high Ha and infrared lu- 
minosities and the anomalously large number of 
Wolf-Rayet (WR) stars in IC 10 (the WR space 
density is highest among the dwarf galaxies, com- 
parable to that in massive spiral galaxies (Massey 
et al. 1992, 2007; Richer et al. 2001; Massey 
and Holmes 2002; Crowther et al. 2003; Vacca 
et al. 2007)) are indicative of a short current star- 
burst affecting the bulk of the galaxy. 
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Lozinskaya and Moiseev (2007) explained the 
formation of the unique synchrotron superbub- 
ble (SS), which was previously associated with 
multiple supernova explosions (Yang and Skill- 
man 1993; Bullejos and Rozado 2002; Rozado 
et al. 2002; Thurow and Wilcots 2005), by a hy- 
pernova explosion. 

In this paper we investigate emission spectra 
of the ionized gas in the regions of current star 
formation in IC 10 based on observations con- 
ducted with the 6-m Special Astrophysical Obser- 
vatory (SAO) telescope and SCORPIO focal re- 
ducer operating in the mode of a long-slit spec- 
trograph. The main goal of our observations is to 
estimate the metallicity for a large number of emis- 
sion regions. Previously, Lequeux et al. (1979), 
Garnett (1990), Richer et al. (2001), and Lee 
et al. (2003) determined the metallicity only in 
three brightest HII regions of the galaxy. 

In the succeeding sections, we describe the ob- 
serving and data reduction techniques, present 
and discuss the results obtained, and, in conclu- 
sion, summarize our main conclusions. 

OBSERVATIONS AND DATA REDUC- 
TION 

The observations of IC 10 were performed with 
the SCORPIO instrument (Afanasiev and Moi- 
seev 2005) operating in the mode of a long-slit 
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Table 1: Summary of observations 



Slit (PA) 


Date 


AA, A 


8\, A 


Texp, s 


Seeing, arcsec 


PAO 


Aug. 17/18, 2006 


3620-5370 


5 


6000 


2.5 


PAO 


Sep. 2/3, 2008 


5650-7340 


5 


3600 


2.1 


PA132 


Aug. 17/18, 2006 


3620-5370 


5 


4800 


1.4 


PA132 


Aug. 18/19, 2006 


5650-7340 


5 


3600 


1.4 


PA268 


Feb. 14/15, 2007 


6060-7060 


2.5 


1200 


1.4 


PA268 


Jan. 15/16, 2008 


3620-5370 


5 


4800 


1.5 


PA45 


Oct. 29/30, 2008 


3650-7350 


9 


3600 


1.0 


PA331 


Oct. 28/29, 2008 


3650-7350 


9 


6000 


0.9 



spectrograph with a slit about 6' in length and 1" 
in width. The scale along the slit was 0.36" 
per pixel. 

For five slit positions, we obtained spectra with 
a resolution from 2.5 to 10 A. Below, the spec- 
trograms are designated as PAO, PA45, PA132, 
PA268, and PA331, in accordance with their posi- 
tion angles. The corresponding spectral range AA, 
spectral resolution <5A, total exposure time Texp, 
and average seeing are given in Table 1. 

The data were reduced in a standard way; the 
spectra of the stars BD+25d4655, GRW+70d5824, 
and GD 71 observed immediately after the object 
at the same zenith distance were used to calibrate 
the energy scale. 

In the spectra being analyzed, we determined 
the continuum level from the underlying stellar 
population based on spline fitting. The inte- 
grated emission-line flux was measured by means 
of single-component Gaussian fitting. To increase 
the signal-to-noise ratio for faint emission regions, 
we performed an averaging over individual HII re- 
gions when processing the spectrograms, with the 
domain of integration being from 2" to 20" in size. 
The errors in the line fluxes were estimated from 
the synthetic spectra that simulated observations 
with the required signal-to-noise ratio. The ranges 
of errors given below in tables and figures corre- 
spond to 3a. 

RESULTS OF OBSERVATIONS 
The Ionized-Gas Emission Spectrum 

The positions of five slits passing through the 
synchrotron superbubble, the bright region of cur- 
rent star formation, and fainter shell structures 



are shown in Figs, la and lb. The spectroscopi- 
cally confirmed WR stars from the lists by Royer 
et al. (2001) and Massey and Holmes (2002) are 
also marked there. For the convenience of iden- 
tification, the slit spectrograms in the figures are 
marked up in arcsec. 

The measurements performed at the points of 
intersection between spectrograph slits allowed the 
actual accuracy of our observations to be esti- 
mated. Table 2 lists the relative line fluxes mea- 
sured from a pair of slit spectrograms at their in- 
tersections. As we see, within the error limits, the 
agreement is largely good. The differences exceed- 
ing the formal observational errors may be related 
to different domains of integration of the fluxes for 
two slit positions on the nebula with a fine struc- 
ture. 

Results of the spectra reduction are pre- 
sented in Fig. 2 and Table 3. Fig. 2 shows 
the distributions of the relative line intensi- 
ties 7([OII]A3727 + 3729A) / 7([OIII]A5007A), 
J([SII]A6717 + 673lA) / /(Ha), 7([NII]A6583A) 
/ /(Ha), and 7([OIII]A5007A) / 7(H/3) along 
slit PA132. The slit is marked up in arcsec along 
the horizontal axis in the lower panel, the numbers 
at the top denote the HII regions from the catalog 
by Hodge and Lee (1990). "M10" designates the 
faint region around the WR star M10 unrecorded 
in the catalog; SS designates the synchrotron su- 
perbubble. As it is seen in Fig. 1, the slit PA331 
(positions 4-11) passes over the faint region near 
the HLllla in the southeast and positions 108- 
115 correspond to the northern periphery of the 
region HL50. To avoid overloading the figures and 
tables, these two regions on slit PA331 are des- 
ignated in Tables 3 and 4 and all figures simply 
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Fig. 1. — Positions of our five long slits on the Ha image of IC 10: (a) the entire galaxy and (b) the bright 
region of current star formation. The long-slit spectrograms are marked up in arcsec. The circles designate 
the spectroscopically confirmed WR stars from Royer et al. (2001) and Massey and Holmes (2002). 
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Fig. 2.— Distributions of the line flux ratios 7([OII]A3727 + 3729A)/J([OIII]A5007A), 7([SII]A6717 + 
6731 A) //(Ha), 7([NII]A6583A)//(Ha), and 7([OIII]A5007A)/7(H/3) along slit PA132. Positions along the 
slits are given in arcsec at the horizontal axis in the lower panel. The numbers at the top of the lower panel 
denote the HII regions from the catalog by Hodge and Lee (1990). "M10" marks the faint region around the 
WR star M10, and "SS" designates the synchrotron supcrbubble. 



as HLllla and HL50, respectively. The slit PA268 
(positions 68-71) passes between the bright neb- 
ulae HL46 and HL48; this region is designated 
below as HL46_48. 

The columns in Table 3 show the designation 
of HII regions from the catalog by Hodge and 
Lee (1990), our long-slit spectrogram utilized for 
the estimations and the used positions on it, the 
relative intensities of the above-mentioned emis- 
sion lines, and the [SIIJA6717/6731 A doublet line 
intensity ratios. 



The interstellar extinction in Table 3 and Fig. 2 
was taken into account. We determined the color 
excess E(B — V) individually for each HII region, 
where possible. The most reliable results were 
obtained from the low-spectral-resolution spectro- 
grams PA45 and PA331. For these spectrograms, 
we made our estimates by comparing the observed 
Ha:H/3:H7 flux ratio with the theoretical one: 
Ha:H/3:i7 7 = 2.86 : 1.00 : 0.47, which, accord- 
ing to Aller (1984), is valid for a typical electron 
density in a star-forming region of 30-300 cm -3 
when the electron temperature is accepted to be 
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Fig. 3.— Diagnostic diagrams of the line flux ratios: 7([OIII]A5007A)/7(H/?) versus 7([SII]A6717 + 
673lA)/7(Ha), versus 7([NII]A6583A)/7(Ha), versus 7( [Oil] A3727 + 3729 A)/7([OIII]A5007A), and versus 
7([NII]A6583A)/7([SII]A6717 + 673lA). The different symbols indicated at the bottom give the relative in- 
tensities averaged over individual Hll regions from the catalog by Hodge and Lee (1990), over the region 
around WR M10, and over the synchrotron superbubble (SS). The uncertainties that exceed 0.1 dex are 
shown as the error bars. 
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Fig. 4. — Diagnostic diagrams with bright > 1.0 • 10 erg s cm , black asterisks) and faint 

(-Pfja — ' er S s_1 cm" 2 , gray diamonds) regions. The measurements in 5-pixel (l"-8) fields along 

the slits PA0, PA132, and PA268 are presented. The heavy cross and diamond mark the galaxy-averaged 
values found by Hidalgo-Gamez (2005) for bright HII regions and the diffuse ionized gas (DIG). 



11000 K (see discussion below). For the regions 
that did not fall in these slits, we assumed either 
the mean value of E{B — V) = 0.95 mag or the 
color excess for the close region that fell into the 
slits PA 45 or PA 331. Our estimates of E(B - V) 
for different regions in IC 10 range between 0.8 
and 1.1 mag. 

Diagnostic Diagrams of Relative Line Inten- 
sities 

Fig. 3 shows the diagnostic diagrams of rel- 
ative line intensities for IC 10 that are tradi- 
tionally used to compare observations with com- 
puted ionization models: 7([OIII]A5007A) / J(H/3) 
versus J([SII]A6717 + 6731 A) / /(Ha), versus 
7([NII]A6583l) / /(Ha), versus 7([OII]A3727 + 
3729A) / 7([OIII]A5007A), and versus J([NII]A6583A) 
/ 7([SII]A6717 + 673lA). Different symbols indi- 
cate the relative intensities averaged over indi- 
vidual HII regions from the catalog by Hodge 
and Lee (1990), over the region around WR M10, 
and over the synchrotron superbubble (SS). The 
observing uncertainties greater than 0.1 dex are 
shown as the error bars in Fig. 3. 

In Figs. 4a and 4b, the data were combined sep- 
arately by bright (-Fjja — 1-OxlO -15 ergs" 1 cm" 2 ) 
and faint (Fjja - 1-0 x 1CT 15 erg s" 1 cm" 2 ) emis- 



sion regions. The measurements along the 1.8" 
wide slits PAO, PA132, and PA268 are shown in 
Fig. 4. The crosses and diamonds in these figures 
also indicate the galaxy-averaged values found by 
Hidalgo-Gamez (2005) for bright HII regions and 
faint regions of diffuse ionized gas (DIG) (see the 
"Discussion" section). 

Estimates of the Relative Oxygen, Nitro- 
gen, and Sulfur Abundances 

The spectra enable us to estimate the relative 
oxygen abundance in a large number of galactic 
HII regions, both bright and faint. Previously, 
such estimates were made for only three brightest 
regions: HL111, HL45, and the compact clump 
HL106a inside HL106 (Lequeux et al. 1979; Gar- 
nctt 1990; Richer et al. 2001; Lee ct al. 2003). 

For the abundance analysis, we employ equa- 
tion (24) from Pilyugin and Thuan (2005), which 
is the revisited calibration dependence of the oxy- 
gen abundance 12 + log(0/H) on the flux ratios 
of strongest [Oil] and [OIII] lines. The classi- 
cal method of estimation from the so-called oxy- 
gen abundance indicator R23 suggested by Pagel 
et al. (1979) and widely used previously was mod- 
ified in a series of papers by Pilyugin (see Pilyugin 
and Thuan (2005)) and is based on two parame- 
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Fig. 5. — Distributions of abundances 12 + log(0/H) (a), log(N + /0) (b), and log(S + /0) (c) along our five 
slits. The HII region identifications are signed at the horizontal axis of the figures. The oxygen abundances 
12 + log(0/H) found using the methods by Izotov et al. (2006), Pilyugin and Thuan (2005), and Chariot 
and Longhetti (2001) are shown by dots with the solid, dotted, and dashed error bars, respectively. 



tcrs: i?23 and the excitation parameter P. Forver- Izotov et al.(2006) and by Chariot and Longhetti 

ification purposes, we also incorporate other tech- (2001). 

niques for the abundance estimation developed by The distributions of abundances 12 + log(0/H) 
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Fig. 5.— (Contd.) 

along the five long slits obtained are shown in 
Fig. 5a and presented in Table 4. The mean 
oxygen abundances in the bright HII regions are 
shown by dots with the error bars. 

As was noted by Pilyugin and Thuan (2005) 
their (24) relation is most reliable for P in the 



range from ~ 0.55 to ~ 1. Having eliminated the 
points corresponding to P < 0.55, we verified that 
this did not change significantly our averaged es- 
timates of the oxygen abundance but only slightly 
increased the measurement accuracy. 

To estimate the relative nitrogen and sul- 
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Fig. 5.— (Contd.) 

fur abundances in various HII regions, we used 
Eqs. (6) and (8) from Izotov et al. (2006) 

The resulting sulfur and nitrogen abundances in 
our regions of IC 10 are shown in Fig. 5a and 5b, re- 
spectively. Note that we actually estimate the S + 
abundance, because we disregarded the S 2+ emis- 



sion in the [SIIIJA6312 line that we can measure 
only in HLlllc. This is also true for the nitrogen 
abundance. 

To properly estimate the N/O and S/O ra- 
tios, we also found the relative oxygen abundance 
12 + log(0/H) using Eqs. (3) and (5) from Izotov 
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Table 2: Emission line flux ratios at the points of intersection of our slit spectrograms 



Spectrogram log(/([OIII] A5007 A) /J(H/?)) log(/([SH] A6717 + 6731 A) /I(tta)) 





1. 


Intersection of PAO, PA132, PA331 




PAO 




0.32 ±0.03 


-0.86 ±0.02 


PA132 




0.30 ±0.03 


-0.68 ±0.04 


PA331 




0.31 ±0.03 


-0.84 ±0.03 




2. 


Intersection of PAO, PA268 




PAO 




0.55 ±0.03 


-1.10 ±0.04 


PA268 




0.55 ±0.02 


-1.25 ±0.03 




3. 


Intersection of PA268, PA331 




PA268 




0.22 ±0.04 


-0.45 ±0.03 


PA331 




0.18 ±0.05 


-0.47 ±0.04 




4. 


Intersection PA132, PA268 




PA132 




0.02 ±0.04 


-0.18 ±0.08 


PA268 




0.13 ±0.05 


-0.33 ±0.04 




5. 


Intersection of PA45, PA331 




PA45 




0.44 ±0.11 


-0.79 ±0.04 


PA331 




0.48 ±0.04 


-0.85 ±0.03 




6. 


Intersection of PA45, PA132 




PA45 




0.61 ±0.21 


-1.06 ±0.11 


PA132 




0.48 ±0.07 


-0.81 ±0.23 



ct al. (2006) 

The characteristic electron density in galactic 
ionized-gas regions that we found from the [SII] 
line flux ratio varies within the range from 30 to 
300-400 cm~ 3 and these variations have virtually 
no effect on the abundance estimates. 

We assume the electron temperature to be 
T e (OIII) = 11000 K in all HII regions. There 
are three reasons for choosing this value. First, 
it is close to a temperature of about 11600 K 
that we determined from the oxygen line flux ra- 
tio in HLlllc. This is the only bright HII region 
where we can detect the [OIII]A4363 A auroral 
line. Unfortunately, the line flux has a large un- 
certainty, because the line is weak. Second, the 
adopted value lies between the ones obtained by 
Lcqucux ct al. (2003) and Lee et al. (2003) for the 
three brightest HII regions in the galaxy. Third, 
this electron temperature provides the best agree- 
ment between the oxygen abundances derived ac- 
cording to Pilyugin and Thuan (2005), and Izotov 
et al.(2006) The temperature in the [Oil] emission 
region was found from the relation (Izotov et al. 
2006) 



t(OII) = -0.744 + t(OIII) • (2.338 - 0.610t(OIII)). 

We assume the same temperature for the [SII] 
and [Nil] emission region. 

The obtained results are collected in Table 4. 
Its columns show the designation of the HII region 
according to the catalog by Hodge and Lee (1990), 
the spectrogram used for our estimates and the 
positions on it, and the derived mean oxygen, ni- 
trogen, and sulfur abundances. The oxygen abun- 
dances estimated using Eq. (24) from Pilyugin 
and Thuan (2005) are denoted in Table 4 by (1). 
Those estimated with Eqs. (3) and (5) from Izotov 
et al. (2006) are denoted by (2). 

We did not average the abundances found from 
different spectra for the same HII region, since 
they correspond to different slit positions on the 
nebula and could differ due to different physical 
conditions: the density and porosity of interstel- 
lar medium, loci of the radiation sources, etc. 

The lower rows in Table 4 show the abundances 
in two regions: HL 111 and HL 45 which cor- 
respond to the regions 1 and 2, from Fig. 1 by 
Lequeux et al. (1979). We derived the abun- 
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Table 3: Emission line flux ratios 



Region 


Sht/position 


/([(Jlll]5007) 
/(H/3) 


/([Sll]6717+6731) 
/(Ha) 


/([IN 11)6583) 
/(Ha) 


/([Sll]6717 
/([S11J6731 


/([(Jll]3727+3729) 
/([(Jlll]5007) 


HLllla 


PAO [-13, -9] 


1.90 ±0.03 


0.178 ±0.002 


0.134 ±0.003 


1.39 ±0.03 


2.64 ±0.33 




PA331 [4, 11] 


2.09 ±0.12 


0.201 ±0.003 


0.099 ± 0.003 


1.10 ±0.01 


1.44 ±0.22 


HLlllb 


PAO [-20, -15] 


1.51 ±0.05 


0.183 ±0.03 


0.148 ±0.003 


1.47 ±0.05 


3.63 ±0.41 


HLlllc 


PAO [-4, 3] 


3.52 ±0.02 


0.081 ±0.001 


0.069 ±0.001 


1.41 ±0.03 


0.80 ±0.10 




PA268 [-6, 5] 


3.47 ±0.01 


0.063 ±0.001 


0.047 ±0.001 


1.40 ±0.01 


0.50 ±0.02 


HLllld 


PA268 [-12,-8] 


2.95 ±0.01 


0.137 ±0.001 


0.097 ±0.001 


1.41 ±0.02 


1.16 ±0.08 


HLllle 


PAO [-7, -5] 


2.44 ±0.06 


0.138 ±0.002 


0.101 ±0.002 


1.50 ±0.05 


1.62 ±0.41 


HL106 


PAO [-42, -21] 


2.89 ±0.07 


0.161 ±0.004 


0.117 ±0.004 


1.46 ±0.08 


1.40 ±0.41 




PA132 [62, 80] 


2.31 ±0.03 


0.191 ±0.007 


0.102 ±0.006 


1.26 ±0.09 


1.23 ±0.08 




PA331 [-20, -8] 


2.42 ±0.11 


0.144 ±0.002 


0.066 ±0.02 


1.02 ±0.01 


1.11 ±0.41 


SS 


PAO [-60, -55] 


0.75 ±0.10 


0.713 ±0.017 


0.247 ±0.013 


1.58 ±0.08 


4.24 ±0.44 




PA132 [84, 100] 


1.33 ±0.13 


0.745 ±0.100 


0.230 ±0.056 


1.24 ±0.26 


3.49 ±0.42 




PA331 [-40, -27] 


1.87 ±0.10 


0.625 ±0.006 


0.207 ±0.04 


1.63 ±0.09 


1.99 ±0.43 


HL37 


PA45 [-12, -6] 


3.34 ±0.24 


0.073 ±0.002 


0.033 ±0.012 


1.55 ±0.05 


0.87 ±0.26 




PA132 [-52, -44] 


2.85 ±0.05 


0.180 ±0.030 


0.064 ± 0.024 


0.99 ±0.22 


0.89 ±0.13 


HL45 


PA45 [-4, 3] 


4.86 ±0.09 


0.079 ± 0.005 


0.045 ± 0.002 


1.12 ±0.03 


0.22 ±0.05 


HL50 


PA45 [12, 27] 


3.44 ±0.08 


0.073 ±0.005 


0.045 ± 0.002 


1.35 ±0.05 


0.69 ±0.14 




PA331 [108, 115] 


3.05 ±0.08 


0.162 ±0.009 


0.086 ± 0.002 


1.49 ±0.14 


1.11 ±0.10 


HL100 


PA132 [49, 58] 


1.65 ±0.02 


0.364 ±0.020 


0.179 ±0.023 


1.24 ±0.11 


2.69 ±0.23 


HL89 


PA132 [35, 43] 


1.32 ±0.2 


0.446 ± 0.029 


0.168 ±0.039 


1.19 ± 0.12 


4.39 ±0.32 


M10 


PA132 [20, 30] 


1.49 ±0.06 


0.753 ±0.065 


0.165 ±0.026 


1.03 ±0.13 


6.41 ±0.53 


HL67 


PA132 [5, 15] 


1.61 ±0.06 


0.578 ±0.076 


0.142 ±0.044 


0.89 ±0.15 


4.13 ±0.48 


HL41 


PA132 [-43, -39] 


1.58 ±0.05 


0.254 ±0.011 


0.126 ±0.013 


1.07 ±0.09 


0.43 ±0.41 


HL36 


PA132 [-61, -55] 


2.86 ±0.03 


0.137 ± 0.019 


0.065 ± 0.007 


0.94 ±0.22 


0.60 ±0.17 


HL22 


PA132 [-116, -108] 


1.53 ±0.03 


0.158 ±0.011 


0.126 ±0.007 


1.17±0.13 


2.60 ±0.37 


HL46.48 


PA268 [68, 71] 


1.49 ±0.08 


0.320 ±0.011 


0.200 ±0.010 


1.66 ±0.09 


2.94 ±0.43 



dances for the relative line intensities measured by 
Lequeux et al. (1979). In both cases the slit in this 
paper was in the east-west direction. The region 
over which the relative abundances were averaged 
spans 3". 8 x 12". 4. Since no localization of the 
spectrograms is given by Lequeux et al.(1979), we 
can make the comparison only approximately. 

DISCUSSION 
Interstellar Extinction 

Variations of the color excess that we revealed 
here for different galactic regions lie within the 
range 0.8 - 1.1 mag. This is considerably smaller 
than the difference from literature, where it ranges 
from E(B - V) = 0.47 mag to E(B - V) = 2.0 
mag (see Sakai et al. 1999; Demers et al. 2004; 
and references therein). In the three brightest 
HII regions, Lee et al. (2003) found E(B - V) = 



0.79 - 1.2 mag, which is similar to our estimates. 
Variations of E(B — V) that we observed can be 
naturally explained by two factors. First, a lo- 
cal higher extinction is possible, since the bright- 
est star-forming region is projected to the dens- 
est cloud of neutral hydrogen and molecular gas 
in the galaxy, and hence can be partly embedded 
into it. A similar effect was also observed by other 
authors. In particular, the extinction estimated 
by Borissova et al. (2000) from the flux ratio of 
the Br7 and Ha lines in the two brightest nebu- 
lae HL111 and HL45 is considerably higher than 
that found from stars. Second, a local lower ex- 
tinction is possible if the stellar wind and SN ex- 
plosions in young groups sweep out the surround- 
ing dense gas. In particular, in the central star- 
forming region Vacca et al. (2007) found E(B — V) 
— 0.6 mag from massive blue stars and concluded 
that the young stars born in the last star forma- 
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Table 4: Relative O, N + , and S + abundances 



Region 


Slit/position 


12 + log(0/H) (1) 


12 + log(0/H) (2) 12 + log(N+/H) 


12 + log(S + /H) 


HLllla 


PAO [-13, -9] 


8.34 ±0.24 


8.28 ±0.12 


6.26 ±0.13 


5.46 ±0.12 




PA331 [4, 11] 


8.05 ±0.34 


8.15 ±0.15 


6.54 ±0.16 


6.00 ±0.16 


HLlllb 


PAO [-20, -15] 


8.27 ±0.22 


8.29 ±0.13 


6.38 ±0.13 


5.55 ±0.13 


HLlllc 


PAO [-4, 3] 


8.24 ±0.15 


8.24 ±0.10 


5.92 ±0.10 


5.07 ±0.10 




PA268 [-6, 5] 


7.94 ±0.10 


8.15 ±0.10 


6.05 ±0.12 


6.56 ±0.22 


HLllld 


PA268 [-12,-8] 


8.28 ±0.11 


8.05 ±0.16 


6.10 ±0.18 


5.33 ±0.19 


HLllle 


PAO [4-7, -5] 


8.29 ±0.38 


8.24 ±0.18 


6.21 ±0.19 


5.42 ±0.18 


HL106 


PAO [-42, -21] 


8.41 ±0.42 


8.28 ±0.22 


6.21 ±0.22 


5.41 ±0.22 




PA132 [62, 80] 


8.06 ±0.14 


8.15 ±0.10 


6.01 ±0.10 


5.32 ±0.10 




PA331 [-20, -8] 


8.03 ±0.42 


8.14 ±0.27 


6.39 ±0.27 


5.88 ±0.28 


SS 


PAO [-60, -55] 


7.59 ±0.42 


8.04 ± 0.43 


6.44 ±0.41 


5.97 ±0.43 




PA132 [84, 100] 


8.07 ±0.18 


8.22 ±0.21 


6.45 ±0.38 


6.00 ±0.41 




PA331 [-40, -27] 


8.12 ±0.42 


8.19 ±0.39 


6.88 ±0.39 


6.50 ±0.39 


HL37 


PA45 [-12, -6] 


8.23 ±0.42 


8.23 ±0.20 


6.23 ±0.24 


5.63 ±0.21 




PA132 [-52, -44] 


8.07 ±0.23 


8.17 ±0.10 


5.96 ±0.22 


5.39 ±0.17 


HL45 


PA45 [-4, 3] 


7.92 ±0.17 


8.21 ±0.10 


6.24 ±0.12 


5.60 ±0.10 


HL50 


PA45 [12, 27] 


8.10 ±0.36 


8.20 ±0.13 


6.32 ±0.14 


5.58 ±0.13 




PA331 [108, 115] 


8.30 ±0.16 


8.25 ±0.10 


6.56 ±0.10 


5.80 ±0.10 


HL100 


PA132 [49, 58] 


8.16 ±0.16 


8.23 ±0.10 


6.31 ±0.15 


5.65 ±0.10 


HL89 


PA132 [35, 43] 


8.23 ±0.13 


8.29 ±0.10 


6.30 ±0.27 


5.72 ±0.10 


M10 


PA132 [20, 30] 


8.74 ±0.36 


8.49 ±0.19 


6.43 ±0.29 


6.04 ±0.20 


HL67 


PA132 [5, 15] 


8.52 ±0.23 


8.36 ±0.13 


6.34 ±0.20 


5.86 ±0.17 


HL41 


PA132 [-43, -39] 


7.73 ± 0.42 


7.79 ±0.39 


6.24 ±0.39 


5.55 ±0.39 


HL36 


PA132 [-61, -55] 


7.86 ± 0.42 


8.09 ±0.15 


5.96 ±0.19 


5.28 ±0.19 


HL22 


PA132 [-116, -108] 


8.28 ±0.12 


8.18 ±0.14 


6.25 ±0.16 


5.44 ±0.16 


HL46.48 


PA268 [68, 71] 


8.10 ±0.40 


8.21 ±0.21 


6.58 ±0.23 


5.85 ±0.21 






from observations by Lequeux et al. (1979) 






HL111 




7.74 


8.17 


6.23 


5.48 


HL45 




7.92 


8.31 


6.24 


5.40 



tion episode are located at the front boundary of 
the above-mentioned dense cloud in front of the 
old stellar population. 

In this paper we present one of the most de- 
tailed studies of the oxygen, nitrogen, and sul- 
fur abundances in IC 10. Previously, Lequeux 
et al. (1979) found 12+log(0/H) = 8.17 for HL111 
and 12 + log(0/H) = 8.45 for HL45. The lower 
rows in Table 4 give our estimates based on the 
measurements by Lequeux et al. (1979). Accord- 
ing to Garnett (1990), 12 + log(0/H) = 8.19 
for HL45. According to Richer et al. (2001), 
12 + lg(0/P±) = 8.23 ±0.09 in the brightest nebula 
HLlllc. For the brightest objects with detected 
[OIILU4363 A line, Lee et al. (2003) found that 
12 + log(0/H) = 8.32±0.14 in HLlllc, 8.05±0.18 



in HLllle, 12 + log(0/H) > 7.96 in HLlllb, and 
12 + log(0/H) > 7.60 in the clump HL106a in- 
side HL106. In the latter two objects, only the 
upper abundance limit was estimated. As we see, 
all listed results are in a good agreement with our 
observations of HII regions in IC 10. 

As our measurements showed, the relative oxy- 
gen line intensities and, accordingly, the relative 
oxygen abundances vary in a wide range, from 7.6 
to 8.5. To understand whether it is related to the 
observational errors or to the actual differences be- 
tween HII regions, we presented the relative inten- 
sities found, respectively, for bright and faint ob- 
jects in Figs. 4a and 4b. As follows from the figure, 
the measurement errors are slightly larger for the 
faint regions than that in the bright ones. How- 
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ever, the scatter of values in the faint regions does 
not exceed significantly that for the bright ones. 
Therefore, we conclude that the oxygen abundance 
variations in the galaxy are most likely real. At 
the same time, the systematic difference between 
the relative intensities in the bright and faint re- 
gions is observed (see below). 

The relative oxygen abundance averaged over 
all galactic regions that wc investigated here is 
12 + log(0/H) = 8.17 ± 0.35, which corresponds 
to the metallicity Z = 0.18 ± O.14Z . 

Diagnostic Diagrams of Relative Line Inten- 
sities 

The relative intensities of diagnostic lines in the 
bright nebulae and faint regions of diffuse ion- 
ized gas (DIG) in the whole galaxy and over its 
four fields taken apart were estimated previously 
by Hidalgo-Gamez (2005) from nine spectrograms 
combined into one wide band. According to this 
paper, the galaxy-averaged relative intensities are 
/([O///]5007A)//(H/?) = 1.86, /([NII]A6583A) / 
/(Ha) = 0.16, and 7([SII]A6717A) / /(Ha) = 0.28. 
The values averaged over the bright nebulae are 
/([OIII]A5007A) / /(H/3) = 2.05, Z([NII]A6583A) 
/ /(Ha) = 0.12, /([SIIJA6717A) / /(Ha) = 0.19. 
The values averaged over the faint diffuse ionized 
gas (DIG) are /([OIII]A5007A) / /(H/3) = 1.66, 
/([NIIJA6583A) / /(Ha) = 0.19, /([SII]A6717A) / 
/(Ha) = 0.31. Our measurements of these line 
ratios in individual HII regions are in good agree- 
ment with the values above, given the observa- 
tional errors. 

As follows from Figs. 4a and 4b, similar to other 
Irr galaxies investigated in detail, the diagnostic 
diagrams slightly differ between the bright and 
faint regions of IC 10: the relative intensities of 
[SII] lines are appreciably higher in faint regions as 
well as those of [Nil] lines are also slightly higher, 
whereas /([OIIIJA5007 A)//(H /3 ) weakens in these 
regions. The similar picture is also observed in 
spiral galaxies, for example, in M31 (see Galarza 
et al. 1999). The value of /([SII])//(Ha) in the 
bright compact HII regions is lower than that in 
the faint diffuse and ring nebulae, not to mention 
the DIG. The noted differences can be explained 
by a decrease in the ionization parameter when 
transiting from the compact to faint diffuse re- 
gions. This effect is predicted by photoionization 
models. The contribution from the gas emission 



behind the front of the shocks triggered by super- 
nova explosions and stellar wind also leads to the 
same effect of the relative intensity Z([SII])//(Ha) 
enhancement. The filamentary structure of the 
entire Ha emission region in the galaxy suggests 
that the shocks play a prominent role in IC 10. 

It should be noted that the constructed diag- 
nostic lines diagrams agree poorly with the cur- 
rently available photoionization models for the 
metallicity Z = 0.2Z Q found above from strong 
oxygen lines. In particular, the observed diag- 
nostic diagrams are consistent with the families of 
theoretical diagnostic curves (Dopita et al. 2006) 
only for a metallicity from Z = OAZq to Z= (1 - 
2) Z Q . 

In Fig. 6, our data are compared with the re- 
sults of Cid Fernandes et al. (2007) and Asari 
et al. (2007), who summarized the gas metallicity 
estimates for ~ 85000 star burst galaxies. On the 
/([OIII]A5007A)//(///3) versus /([NII]A6583A)//(Ha) 
diagram, these authors give the number of objects 
in each sector of the galaxy localization band and 
the mean metallicity that corresponds to the sec- 
tor. As we see, the HII regions in IC 10 fall 
nicely on the galaxy localization band but fall 
within the range of mctallicities from Z = O.SZq 
to Z = 0.6Zq. One distinctive region in IC 10, the 
synchrotron superbubblc, falls near Z = 0.8Z Q . 
This object is discussed below. Note, however, 
that Asari et al. (2007) recorded a systematic 
shift by 0.2 dex between their metallicity esti- 
mates and those of Pilyugin and Thuan (2005) 
and Izotov et al. (2006) for 177 "common" objects 
made from the relations that we used here. It 
turns out that the lower the metallicity, the larger 
the shift. For this reason, Asari et al. adopted the 
metallicity Z = 0.2Z Q for Sector A rather than 
Z = 0.3 Zq and therefore Z = 0.2 Z & is the lower 
boundary of the range into which our regions of 
IC 10 fail. 

It is also interesting to compare the observed 
diagnostic diagrams with the theoretical photoion- 
ization models by Chariot and Longhetti (2001), 
since these authors separately took into account 
the changes in each of the parameters defining the 
model calculations, such as the cluster age, the 
initial mass function (IMF), ionization parameter, 
density of interstellar medium, and metallicity. 

Comparison of the diagnostic diagrams for IC 10 
with the distributions of /([OIII]//(H/3) as a 
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Fig. 6. — Comparison of our line intensity ratios with results by Cid Fernandes et al. (2007) and Asari 
ct al. (2007). The mean metallicity in each of the gray sectors is shown to the right of the localization 
band of 85000 starburst galaxies. The circle and the square indicate the galaxy- averaged values found by 
Hidalgo-Gamez (2005) for bright HII regions and the diffuse ionized gas (DIG), respectively. 



function of J([OII])/J([OIII]), I([SII])//(Ha), 
and 7([NII])/7([SII]) constructed by Chariot 
and Longhctti (2001) for galaxies of various 
types and HII regions showed an excellent agree- 
ment between them and our observing line flux 
ratios. In Fig. 7 our observations are com- 
pared with the models computed by Chariot and 
Longhctti (2001) for various metallicities and var- 
ious effective ionization parameters. The large 
square, circle, and triangle on our diagnostic di- 
agrams (left panels) mark the regions of reduced 



(Z = 0.2Zq), solar [Z = Z@), and enhanced 
(Z = 2Zq) metallicity from models by Chariot 
and Longhctti (2001), see Figs. 1, 2, and 3 in 
their paper. The large square, circle, and trian- 
gle on the right panels indicate three values of 
the model flux ratios for the effective ionization 
parameter logt/o of —3.0, —2.5, and —2.0 from 
the calculations by Chariot and Longhetti (2001). 
As we see, in complete agreement with the model 
calculations, the J([OIII]A5007A)//(H/3) versus 
I([NII]A6583A)//([SII]A6717 + 673lA) dia gnostic 
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Different Z 



Different U„ 




-1.5 -1.0 -0.5 0.0 

log([SII]X6717 + 6731/Ha) 




-1.5 -1.0 -0.5 0.0 

log([SII]X671 7 + 6731 /Ha) 





-7 



-10 12 

log([OII]X3727 + 3729/[OIII]X5007) 




-1.0 



-0.5 0.0 0.5 

log([NII]X6583/[SII]X67 17 + 6731) 



-10 12 

log([0ll]X3727 + 3729/[OIII]X5007) 




-2 t 



-0.5 0.0 0.5 1.0 

log([NII]X6583/[SII]A671 7 + 6731) 



Fig. 7. — Comparison of the line flux ratios on diagnostic diagrams with the model calculations by Chariot 
and Longhetti (2001). The large squares, circles, and triangles in the left panels (a) mark the regions of 
reduced [Z — 0.2Z©), solar (Z = Zq), and enhanced (Z — 2Zq) metallicity, respectively, in accordance with 
the model calculations. The same symbols in the right panels (b) indicate the model calculations for the 
effective ionization parameters logj/o — 3.0, —2.5, and —2.0, respectively. 
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diagrams are most promising for estimating the 
metallicity of the galactic gaseous medium. Posi- 
tions of our HII regions in IC 10 on the diagnostic 
diagrams correspond to the metallicities within 
the range of Z= (0.2 - 1) Zq. The dependences 
of 7([OIII]A5007A)/7(H/3) on 7([OII])/7([OIII]) 
and of 7([OIII]A5007A)/7(H/3) on 7([SII]A6717 + 
+673lA)/7(Ha) are well explained by variations 
in the ionization parameter. 

The differences between the metallicity of in- 
dividual HII regions in the galaxy found from 
the i?23 parameter and the metallicity determined 
by comparing diagnostic diagrams with theoret- 
ical photoionization models can be explained by 
a number of factors. First of all, the photoioniza- 
tion models are computed for rich star clusters. In 
particular, Dopita et al. (2006) assumed a cluster 
mass of M > 10 3 M Q , the IMF by Miller and Scalo 
ranges from 1 to 100M Q , while the star clusters 
in IC 10 are comparatively poor (see Hunter 2001). 

In addition, the contribution from the gas emis- 
sion behind the front of the shocks can be sig- 
nificant in IC 10, as suggested by the filamen- 
tary structure of the entire star-forming region in 
this galaxy. Actually, the term "diffuse ionized 
gas" (DIG) can be applied to IC 10 only condition- 
ally, because the faint emission regions are char- 
acterized mostly as filamentary shell and/or arc 
structure. Therefore, the faintest regions inside a 
number of shells that appear diffuse most likely 
represent the front or back sides of the shells. 

The synchrotron supcrbubble in IC 10 is a 
unique object and deserves a more detailed discus- 
sion. According to Lozinskaya and Moiseev (2007) 
and Lozinskaya et al. (2008), the image in [SII] 
lines better than in Ha reveals the optical shell 
identified with the synchrotron superbubble. The 
expansion velocity of the bright knots and fila- 
ments in the optical shell found in these papers, 
50-80 km s _1 , confirms the estimate by Ramsey 
et al. (2006). The gas mass is M ~ 4 x 1O 5 M 
and the kinetic energy is E^- m ~ (1 — 3) ■ 10 52 erg. 
This amount of released energy corresponds to the 
explosions of ten supcrnovae plus the stellar wind 
from their host association. Alternatively, the su- 
perbubble may be a result of a hypernova explo- 
sion, as was suggested by Lozinskaya and Moi- 
seev (2007). The kinematic age of the superbub- 
ble, t ~ (3 - 7)-10 5 yr, is a strong argument for a 
hypernova, because it takes more than t ~ 10 7 yr 



for the explosions of ten supernovae in a local 
galactic region. 

As it can be seen in Fig. 1, the regions from 85" 
to 120" on spectrogram PA132, from -60" 
to —25" on spectrogram PA331, and from —50" 
to —80" on spectrogram PA0 (the latter passes 
over the very edge of the superbubble) correspond 
to the synchrotron superbubble. 

Our measurements indicate that the relative 
intensity 7([SII]A6717 + 673lA)/7(Ha) in the 
region of the synchrotron superbubble reaches 
0.7-0.75, which strongly suggests a gas emis- 
sion behind the shock front (Dopita and Suther- 
land 1995, 1996). This is also suggested by 
the ratios of 7([NII]A6583A)/7(Ha) ~ 0.25 and 
7([OII]A3727 + 3729A)//([OIII]A5007A) ~ 2 - 4 
that are higher than those in the nearby bright 
region HL106, and also by the lower relative in- 
tensity 7([OIII]A5007A)/7(H/3). 

The region of the synchrotron superbubble on 
the diagnostic diagrams also differs from HII re- 
gions by a lower ratio 7([OIII]A5007A)/7(H/3) and 
by high 7([SII]A6717 + 673lA)/7(Ha), 
7([NII]A6583A)/7(Ha), and 
7([OII]A3727 + 3729A)/7([OIII]A5007A). This 
suggests the dominance of the collisional exci- 
tation and is typical in old supernova remnants. 

Lozinskaya and Moiseev (2007) found a mean 
n e ~ 20 - 30 cm' 3 from the intensities of the 
[SII](A6717/673lA) lines for the brighter northern 
part of the superbubble (the region 90" - 105" 
on PA132), where the errors are relatively small. 

For the region —50" 60" on spectrogram PA0, 

we find a similar mean density n e ~ 30 - 40 cm~ 3 . 
In the southern faint part of the superbubble, in- 
dividual knots can be distinguished in the spec- 
trogram PA 132 in which the density reaches n e ~ 
200 - 300 cm~ 3 , but the accuracy of these esti- 
mates is low. The gas density northward of the 
synchrotron superbubble (the positions range be- 
tween 80" - 90" on PA132) is n e ~ 200 cm~ 3 . 

The oxygen abundance estimate in the region 
of the synchrotron superbubble is ambiguous. In 
general, the present-day models of the gas emis- 
sion spectrum behind the shock front developed by 
Allen et al. (2008) allow the abundances of heavy 
elements, in particular oxygen, to be determined 
from the available diagnostic relative line intensi- 
ties if the ambient gas density and the shock veloc- 
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ity are known. However, in our case the problem is 
complicated by the presence of the WR star M17 
in the central region of the synchrotron superbub- 
blc, whose ionizing radiation may be significant 
or even dominant in close neighborhoods. On the 
other hand, we cannot properly use the technique 
applied above to HII regions without taking into 
account the gas emission behind the shock front. 

A careful calculation including both emission 
mechanisms is beyond the scope of this paper. Be- 
low we formally estimate the oxygen abundance 
in the synchrotron superbubble independently by 
two methods, realizing that both arc improper. 
Employing the R23 parameter, we found the abun- 
dances 12 + log(0/H) = 8.07 ± 0.18 from spectro- 
gram PA132, 7.59±0.42 from PAO, and 8.12±0.42 
from PA331 using the technique of Pilyugin and 
Thuan (2005) or, respectively, 8.22 ± 0.21 from 
spectrogram PA132, 8.04 ± 0.43 from PAO, and 
8.19 ± 0.39 from PA331 using the technique of Izo- 
tov et al. (2006) (see Table 4). Since the mea- 
surements were made for three different regions of 
a very faint extended object, these independent 
estimates of the oxygen abundance are in a very 
good agreement. 

A detailed kinematic study of the synchrotron 
superbubble (Lozinskaya et al. 2008) clearly re- 
vealed weak Ha wings in the velocity range 
from —200 to —450 km s _1 . These high- velocity 
features in the line correspond to the velocity of 
the shock wave triggered by a hypernova explosion 
of about 100-200 km s" 1 . Comparison with the 
models of gas emission behind the front of a shock 
(Allen et al. 2008) propagating with such veloc- 
ity gives the best agreement at the abundance 
12 + log(0/H) ~ 8.15 with the models of purely 
collisional excitation, or 12+log(0/H) ~ 8.35 with 
the models that include the gas pre-ionization by 
a shock. 

CONCLUSIONS 

We present the most detailed to-date study of 
the O, N + , and S + abundances of the individual 
HII regions and the synchrotron superbubble in 
the dwarf Irr galaxy IC 10 that we estimated from 
five spectrograms taken with the SCORPIO in- 
strument at the 6-m SAO telescope. 

The relative O, N + , and S + abundances in in- 
dividual HII regions of the galaxy IC 10 lie within 



the following ranges: 12 + log(0/H) = 7.59 - 8.52 , 
12 + log(N+/H) = 5.92 - 6.58 , and 12 + log(S+/H) 
= 5.00 - 6.00. The mean metallicity of the galaxy's 
gaseous medium is Z = 0.18±0.14Z©. The metal- 
licity estimates from the comparison of diagnostic 
diagrams with photoionization models are found 
to be less reliable than those from the strong oxy- 
gen line ratios. 

The oxygen abundance estimates in the syn- 
chrotron superbubble from the measured relative 
line intensities must be based on the model of the 
combined action of a shock (the result of a hyper- 
nova explosion) and ionizing radiation (the pres- 
ence of a WR star in the central region) . So far the 
estimates have been made separately within the 
framework of these two models. With the model of 
an HII region, we found that 12 + log(0/H) =7.59 
- 8.34. The comparison with the models of gas 
emission behind the front of a shock propagating 
with the velocity of 100-200 km s -1 gives the best 
agreement at the abundance 12 + log(0/H) ~ 8.15 
with the models of purely collisional excitation or 
12 + log(0/H) ~ 8.35 for the models that include 
the gas pre-ionization by a shock. 
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